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Asymptotics of alternating harmonic series with 

attennation 

Sergey Sadov* 


Abstract 

We find the asymptotics of the series exp(—t/n) as 

t —)■ + 00 . The answer is an oscillating function of t dominated by 
exp(—(27rt)^/^). The intermediate step is to find the asymptotics of 
the two-dimensional Fourier transform F{^) of the function F{x) = 

(1 -f exp(||a;||p)-i as |g|| ^ oo. 
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The problem discussed in this note is an yet another example of a chal¬ 
lenge born of a teaching mishap. I offered it by mistake among a set of 
exercises on the Euler-Maclaurin formula in a course of asymptotic analysis 
at Memorial University of Newfoundland in the Fall 2011. 

Problem. Find asymptotics of the series 

n=l 


as t —)■ 4 - 00 . 

In the interests of those who want to take up the challenge the answer is 
only given at the end of the note. The solution below is long but detailed. 
It is intended to be understood by an asymptotic analysys course student 
who wants to get their hand dirty with the Saddle Point Method. 

*E-mail: serge.sadov@gmail.com 
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Solution. 

Part 1: Derivation of integral representation 
Let us consider a more general series 

OO Yi 

n=l 

Lemma 1 (a) The series (§ converges absolutely if \z\ < 1 for any G M. 

(b) Ifiy>l, then the series converges absolutely and uniformly in the closed 
unit circle \z\ < 1. 

(c) If 0 < u < 1, \z\ = \, z the series converges conditionally. Also, 
for any e G (0,2), the convergence is uniform in the region 

De = {z : |- 2 ;| < 1, Re z < 1 — e} . 

(d) Consequently, if \z\ = 1, z 1, then 

S{z,n,t)= lim S{zp,v,t). (3) 

p^i- 


Proof. Parts (a) and (b) are obvious since the series is dominated by 
\z\'^n~'^ in the case (a) and by the z-independent sum ™ case 

(b). 

To prove (c), let us write = 1 — 6n, where 5n = 0{n~^) (for a fixed 
t). The series ^ convereges absolutely and uniformly in {| 2 ;| < 1} 

(similarly to the series S{z,i' + l,t)). It remains to establish the convergence 
properties (c) for the series 

OO 

n=l 


It is a standard application of Dirichlet’s test: the sequence {n decreases 
and the uniform in bound for the partial sums 


N 


E- 

n=l 


2 2 
-r < - 

1 — Z £ 


holds. 

^The function Li,^(2) defined by this series is called polylogarithm of order v. 
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To prove (d), write 2 : = e*®,where 6 G (0, 27r). Let e > 0 be such that 
e < 1 — cos 6. Then pz G for any p G [0,1]. The series 

y- [pe ) ^-t/n 

71=1 


converges uniformly w.r.to p. Hence termwise passing to the limit as p —?• 1 
is justified, and Q follows. □ 


Remark. In Lemma the parameter t can be any complex number. For 
fixed 2 ; (l^l <1, 2 ; 7 ^ 1) and ly > 0, the function S{z, u, t) is an entire analytic 
function of t. 

It is interesting also to note the relations 


dS{z, V, t) 
Wt 


-S{z,v + l,t), 


dS{z, V + l,t) 
dz 

d‘^S{z, V, t) 
dt dz 


= z ^S{z, V, t), 
—z~^S{z, V, t). 


□ 


Recall the series representing Bessel function of order n > 0: 

(-1)^(u/2)2^+’^ 


Jn{u) — 


fc =0 


k\ r(/c + n + 1) ’ 


( 4 ) 


Lemma 2 For i' > 1 and \z\ < 1, the function S{z,n,t) has the integral 
representation 


POO 

S{z,n,t) = - 

Jo J- 


ze 


x\ 


Ji,-i{2\/tx) dx. 


— ze * \t 

The same representation remains valid if \z\ = 1, z In particular, 

COO 


s{t)=s{-i,i,t)=- r 

Jo 


+ 1 


Jo( 2 \/te) dx. 


( 5 ) 


( 6 ) 


Proof. Suppose \z\ < 1. Since 
T{k + n) 


n 


fc+Z/ 


f 


e-n^ xk+''-Ux, 


3 










we have 


S{z,v,t) 


OO OO 




n=l k=0 

OO OO 

E-"E 

n=l k=0 
OO 


k\ 


i-tY 


k\ T{k + v) Jq 


g-nx ^k+u-1 


E 

k=0 

OO 

E 

k=0 

poo 

Jo 


k\ r{k + u) Jq 

i-t)^ r 

kl r{k + Jy) Jo 1 — -ze 


OO 

'.z-e 

\n=l 
^—x 


E- 


n^-nx ^k+v-l 


ze 


k+u—l 


—X 


dx 


ze 

1 — ze~ 


. dx. 


k=0 


k\T{k + i^) 


In the above transformations, convergence is absolute throughout, so we can 
interchange the order of summations and integration freely. 

Setting u = 2^/tx and comparing the integrand in the RHS with Q, we 
see that 


if- {-txf 

^ k\T{k + u) 


\tJ ^ A:!r(A; + i^) 



Ju—l{u) , 


SO that (j^ follows. 

Denote temporarily the RHS of ([^ by S^{z,u,t). The integral in (j^ 
converges absolutely if |z| < 1, z / 1. (Convergence at x —)> oo takes place 
for any z and convergence at x —)• O"'' is certain if z/ > 1, z ^ [l,oo).) 
Consequently, if |z| = 1, z / 1, 


lim 5* {zp, u, t) = S'* (z, u, t). 

p->i- 

On the other hand. Lemma asserts the same for S{z,i',t). Thus, 
remains valid in the case |z| = 1 , z 7 ^ 1 by continuity. □ 


Corollary. Let 

X = 2Vt and S{t) = S^{X). (7) 

Taking ^/x as the new variable of integration in ([^ , we get 

2r 

S*iX) = - ■/o(Ax) 3 dx. ( 8 ) 

Jo + 1 

is called the Hankel transform of the function —2x/{exp{x^) + 1). 
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Part 2: Transformations of the integral ([^ 

Substituting the integral representation 

2-k 

^iu cos 9 



de, 


into Q we express 5* (A) as a double integral 



( 9 ) 


Experimental study suggests that 5*(A) = 0(e Unfortunately, it is 
not clear how a precise asymptotics can be derived from the representation 
(§, because the edge x = 0 interferes in all attempts: integration by parts 
or rotation of the integration path in the complex x-plane. We will not 
describe such attempts in detail, since the successful solution will be based 
on a different transformation of the integral Q. Note only that if instead 
of the factor 2x we had an even function of x in the integrand of , then it 
would be possible to extend domain of integration w.r.t. x to M, getting rid 
of the edge; the saddle point method would then be applicable in a standard 
manner. 

The representation we need will come from a known general formula of 
Fourier analysis. 

Lemma 3 Suppose f^\fi'r')\rdr < oo. Let r = r{x,y) = 

Then the 2-dimensional Fourier transform of F{x,y) = f{r(x,y)) is 


■OO 


F{(,,r]) = / F{x,y)e ^^^dxdy = In Jo{rp) f{r) r dr, (10) 
JM2 Jq 

where p = 

Proof. Using the polar coordinates (r, 6) in the (x, y) plane and {p, 4>) in the 
(?) v) plane, we get 

poo p2'K 

F{^,7])= / / f{r)e-^’^P^°<‘^-^Uerdr. 

Jo Jo 

Clearly, the change of variable 9 ^ 9 — (p shows that the integral does not 
depend on (f>; thus 
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as stated. 


□ 


Corollary. The integral 
integral 


1 poo poo 

s*W = -- / T 

^ J — OO J — OO 1 


can be written as the 2-dimensional Fourier 

iXx 




dx dy. 


( 11 ) 


Indeed, it suffices to replace r in Q by r and to use (10) with /(r) = 
2 /(e-^+i + l), (^,^) = (-A,0). 


Remark. Expanding (1 -|- as a geometric series and integrating 

separately w.r.t. x and y, one obtains an independent proof of the fact that 
S{t) equals the RHS of (11) when A = 2y/i. 


Part 3: Asymptotic analysis of the integral (11) 

Asymptotics of the integral (11) is determined by the complex singular¬ 
ities of the integrand. As a preparation, let us locate complex zeros in the 
2 ;-plane of the denominator 

Q{z,y) = l + e^"^y\ 

where z = x + iu, and x, u ,y are real. 

The equation Q{z,y) = 0 is equivalent to 


z^ + y^ = i{2k -|- l)7r, k € Z. 


Or, in the real form, 


2,2 2 
x +y = u , 


XU = TT { k + - 


( 12 ) 


Precise asymptotic analysis of the integral (11) is significantly more com¬ 
plicated than singularities-based analysis of Fourier integrals in the one¬ 
dimensional case. As a first, rather easy step, we obtain a rough exponential 
o-estimate. 

Lemma 4 For any a < y^7r/2 we have 

,S,(A) = o(e-“^). (13) 
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Proof. Clearly, for any fixed y > 0 and h > 0 


|1 + e 


{x+ihp+y'^ 


OO 


as X ^ ±oo. Therefore, assuming Q{x + ih, y) 7 ^ 0, Vx G M, we have 


roo ^ r 

j-ooQ{x,y) y_ 


00 ^i\{x+ih) 

-00 Q{x + ih,y) 


+ 27ri Res 


^iXz 

Qiz,y)' 


(14) 


where the sum of residues is taken over all z = x + m such that 0 < u < h 
and Q{z, y) = 0. 

The equation Q{x + iu,y) = 0 does not have solutions with x,y G M, 
0 < n < y/TTf2. Indeed, looking at the system (12), we see that |x| < u and 
|7r/2| > |xm| > vf. 

Therefore, for h < Y^7r/2, the sum in the RHS of (14) is void. Thus 

^iX{x^ih) 


1 poo poo 

s*W = — / - 

^ J— OO j— OO 1 


g(x+2/i)2+y^ 


dxdy = 0{e 


-\h\ 


Given a G (0, 7 r/ 2 ), the estimate (13) follows by taking h G (a, 7 r/ 2 ). 

□ 

We want eventually to nail down an asymptotical term of the exponential 
order 0{e~^^^) precisely. 

The procedure consists of simple steps, but it is rather delicate overall. 
We need to take into consideration some poles of Q{z,y)~^ to obtain a 
nonempty the residue sum in the RHS (14). On the other hand, we must 
avoid stepping on a pole during the double integration (in both x and y). 
Moreover, any estimates we get while keeping y fixed should be uniform or 
explicit enough to justify subsequent integration with respect to y. 

The following lemma will allow us to control the integrand in the double 
integral uniformly. 


Lemma 5 Suppose u > 0 and a closed set ifT C M are such that the equation 
Q{x + iu, y) = 0 does not have solutions x G M, y G RT. Then there exists 
a = a{a, K) > 0 such that 

\Q{x,y)\ > Vx G M, y € K. (15) 

Proof. We will estimate \Q{x,y)\ from below in two cases separately: first 
for large x^ + y^, then in a bounded region of the (x, y) plane. 
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1) Suppose that + \n.2. Then hence 

\Q{x,y)\ > - 1| = - e“') 

p-U^ 2 2 

> _ _ pX^+y 

~ 2 


2) The set 

Vt. = {{x,y)\x^ + y‘^ < + \n2, ye K} 

is a bounded closed subset of hence compact. By assumption, Q{x,y) / 
0 when {x,y) G 0. The function /Q{x,y) is continuous in fl, hence 

bounded: /Q{x,y)\ <C. 

The inequality ( [T^ with 

ll 

2 ’ Cj 



follows. □ 

Next comes the crusial step. We choose a splitting parameter b for 
integration with respect to y and write (0 in the form 


- 7 r 5 *(A)= [ T{y,X)dy+ f T{y,X)dy = h + h, 
J\y\>b J\y\<b 


where 


T{y,X) = 


A\x 


_oo 1 + +y 


dx. 


(16) 


(17) 


We are going to make a complex shift x —)> x + m in the integral 0 . 
However, we will use two different values of a in the two cases, \y\ > b and 

\y\ < b. 

In the case \y\ > 6 , we will take a = oi so as to ensure that no zeros of 
Q{z, y) lie in the strip 0 < Imz < ai. It will be possible to choose ai > Y^ 7 r /2 
and to obtain the estimate Ii = 0{e~^°‘^). 

In the case \y\ < b, we will take a = 02 > tt/2 to ensure that the 
equation Q{z,y) = 0 has exactly two solutions z± = z±{y) (with the same 
imaginary part) in the strip 0 < Imz < oi. The principal term of the 
asymptotics of the integral ( |11[ ) will come from the residue part of the RHS 
in (14). 


The three parameters b, ai and 02 in the outlined program are not deter¬ 
mined rigidly and can be varied as long as we don’t cross certain boundaries. 
Details are clarified in the following two lemmas. 
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Lemma 6 For any a satisfying the inequalities 



there exists a unique b > 0 such that 

(i) the equation Q{x + iu,y) = 0 does not have real solutions {x,u,y) with 
0 < u < a and \y\ > b; 

(ii) for every y G {—b,b), the system 


Q{x + zw, y) = 0, 0 < u < a 


has exactly two real solutions {x,u). They have the same imaginary part 
and opposite real parts, that is, the solutions are of the form x = ±x*(y), 
u = 


Proof. We will prove that b = \/cP — Xq, where xq 
pair X = xo, y = b satisfies the system (12) with u 


7 r/(2a). Note that the 
a, k = t). 


(i) If Q{x + iu,y) = 0 and \y\ > b, < u < a, then < 

of — b'^ = Xq. Hence |xzt| < xqu = 7r/2, and the equation xu = TT{2k + 1) 
cannot hold. 


(ii) Let 7 = 7r/2 or —7r/2. Substituting x = 'y/u into the equation 
+ y^ = , we get 

4 2 2 2 n 

u — y u —7 = 0 . 

Since u is real, there must be > 0, hence 


^ + Vy2 + 472) . 


(18) 


Let u*(y) be the (positive) square root of the RHS. Clearly, u*(y) is an 
increasing function. By definition of the number b (at the beginning of the 
proof), we have u*( 6 ) = a. Hence, |y| < b implies 0 < n* < a, as required. 

If we replace the value I 7 I = 7 r /2 by 7 r(/c + l/2) with A: > 1, then the three 
inequalities > 7 ^, 0 < u < a, and < 37 r /2 < I 7 I, are incompatible. 

The rest is trivial. The value 7 = 7 r /2 yileds the solution u = u*, 
X = X* = 7r/(2M*), while 7 = —tt/2 yileds the solution u = u*, x = —x*. 
□ 


Lemma 7 

ities 


There exist real numbers 6 > 0 and oi, 02 satisfying the inequal- 


vr 

- < ai < 02 < 



(19) 
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with the following properites: 

(i) the equation Q{z, y) = 0 does not have solutions (z, y) with y G M, \y\ > b, 
0 < \mz < ai; 

(ii) for every y G {—b,b), the system 

Q{z, y) = 0, 0 < Imz < 02 


has exactly two solutions z±{y) = ±x*(y) + iu^,{y). 


Proof. Choose ai and 02 satisfynig the inequalities (19) arbitrarily. Take 
any a G ( 01 , 02 ) and determine b as in Lemma Let us check that the 
required properties are in place. 

(i) If |y| > b, then by Lemma(i), the equation Q{z,y) does not have 
solutions with 0 < Imz < a. The constraint Imz < oi is even stronger. 

(i) Let |y| < b. By Lemmaj^ (ii), the equation Q{z,y) has two solutions 
z = satisfying the inequality 0 < Imz < a. Relaxing the constraint 

to 0 < Imz < 02 could potentially bring in extra solutions. For every such 
solution, u = Imz would be defined by (18) with I 7 I = n^k + 1/2), k > I. 
Then > Sj27r, which contradicts the condition n < 02 < ^ 37 r/ 2 . □ 


Lemma 8 Let oi, 02 , b be as in Lemma^ and Ii, I 2 as in (16). Then 

R = 0(e-^“i). 


I 2 = -ni ^ 


rb ^iXz±{y) 

-b z±(y) 


dy + 0(e-^“2), 


( 20 ) 


where z±{y) are the two solutions of the equation Q{z,y) = 0 defined in 
Lemma^ 


Proof. 1) Consider the integral (© with 
where we set h = ai, 


> b. By Lemma [^(i) and (14) 



gA(-ai+ia;) 

Q{x + iai,y) 


dx. 


Therefore, 


\h\ < 



|r(y,A)|dy<e-^“i 



dx dy 

\Q{x + iai,y)\' 
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The conditions of Lemma are met with u = oi, K = (— 00 ,— 6 ] U [ 6 , 00 ). 
The estimate (20) for Ji follows at once. 


2) In (14) we set h = a 2 - By Lemma[^(ii) there are two poles, 2 ; = z±{y), 
that contribute to the sum of residues. The residue at 2 ; is 


r,i\Z 


^iXz 


ryiXz 


Q'^{z,y) 2ze^^+y^ -2z' 

since at a pole Q{z,y) = 0 and = —1. We have identified the sum 


in ( 20 ). 


part of the RHS in (14) with the first part in the RHS of the second estimate 


To show that the integral in the RHS of (14) is 0{e it suffices to 


repeat the same argument as in part 1 of this proof with obvious changes: 
u = 02, K = [—b,b]. □ 

To simplify the remaining calculations, note that 

2 :+ = X* + iu*, 2 ;_ = —X* + iu*. 


so 


Therefore 


Lemma 9 


— m 


■E 


IZ- = iz+. 
„iXz± 


Z± 


= 27r Re 


iz+ 


fb i\z+{y) 


_2l/2^1/4 gi(A^V2+7r/8) 


( 21 ) 


( 22 ) 


izPy) 

as X ^ + 00 . 

Proof. The function Re(i 2 :+(y)) = —u^{y) attains its maximum —p7rf2 at 


y = 0, cf. (18). 


The saddle point method tells us that 

nb ^i\z+{y) giA2+(0) 


27r 


J-b iz+{y) iz+iO) \-iXzpO) 

Recalling that z'f = ni — y'^, evaluate: 


1/2 


iz+{0) = (-1 + i) 
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(since Rez+,Imz+ > 0). 

Expanding Z-\.{y) in powers of we find 


z+{y) = (^1 - + 0(/), 


hence 


and 


Thus 


4 ( 0 ) = 


—m 


-i4(0) = 

\/7T 


j-b ^i\z+{y) g(-l+i)A-472 / 27r^/2 \ 

J_f) iz+{y) yAe^^*/^y 

Simplifying, we obtain 


1/2 


Combining (20), (21), (22), we get 


I, ^ _23/2„5/4 cos(A442 + 7r/8) 

2 ^ /I e ^^^2 ’ 


and, using (16) 


S.(A) ~ 

A ^/2 

The final answer in the original notation, cf. 0 , is 


Sit) ~ "/**). 


□ 


(23) 


(24) 
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Plots —e'''V^^S'*(A) vs A 



Remark. From our analysis it is easy to see that the (unwritten) error term 

O (e-^V^A-3/2) . 


in (23) is 


The form in which the answer is presented in (23) or (24), is slightly inaC' 


curate: the error term needs not be dominated by the main term everywhere, 
since the latter becomes zero for some values of A (or t). 

The meaning of the ~ sign is that the error term is smaller by its order of 
magnitude, which is characterized by the non-os dilating (amplitude) factor 
in (23) or (24). Regarding the usage of the ~ sign in this and similar 


situations, I disagree with de Bruijn who expressed negative opinion of such 
usage {Asymptotic methods in analysis, end of Sect. 5.11). 
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